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Abstract 

We address continuous weak linear quantum measurement and argue that it is best understood 
in terms of statistics of the outcomes of the linear detectors measuring a quantum system, for 
example, a qubit. We mostly concentrate on a setup consisting of a qubit and three independent 
detectors that simultaneously monitor three non-commuting operator variables, those correspond- 
ing to three pseudo-spin components of the qubit. We address the joint probability distribution 
of the detector outcomes and the qubit variables. When analyzing the distribution in the limit of 
big values of the outcomes, we reveal a high degree of correspondence between the three outcomes 
and three components of the qubit pseudo-spin after the measurement. This enables a high-fidelity 
monitoring of all three components. We discuss the relation between the monitoring described and 
the algorithms of quantum information theory that use the results of the partial measurement. 

We develop a proper formalism to evaluate the statistics of continuous weak linear measurement. 
The formalism is based on Feynman- Vernon approach, roots in the theory of full counting statistics, 
and boils down to a Bloch-Redfield equation augmented with counting fields. 

PACS numbers: 03.65.Ta, 03.67.Lx, 74.50. +r 
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I. INTRODUCTION 



The theory of quantum measurement, being a foundation of quantum physics, is attract- 
ing more and more attention 1 . Intrinsic paradoxes^ are definitely a main reason for studying 
quantum measurements. More motivation comes from the practical needs to understand 
the real solid-state based devices 3 ' 4 developed for quantum computing^. Measurements in 
solid-state setups may provide access to extra variables that facilitate the read out of the 
quantum information stored in the elementary two-level quantum systems (qubits). The 
concept of continuous weak linear measurement (CWLM), where the interaction between 
the detector and the measured system is explicit and sufficiently weak, has been recently 
elaborated in context of the solid state quantum computing&£>&2ii£. CWLM provides a 
universal description of the measurement process and is based on general linear response 
theory^- . It applies to a large class of linear detectors: From common amplifiers to more ex- 
otic on-chip detectors such as quantum point contact^, superconducting SET transistors^, 
generic mesoscopic conductors^ 4 -, fluxons in a Josephson transmission line to measure a flux 
qubits. 

It is an important feature of CWLM that the (quantum) information is transferred from 
a quantum system being measured — a qubit — to other degrees of freedom: those of the 
detector. The outcome of the measurement is thus represented by the detector degrees of 
freedom rather than those of the qubit. We will address both the statistics of the outcomes 
and joint statistics of the outcomes and the qubit degrees of freedom. 

We stress the difference between the detector outcomes and the outcomes of a projective 
measurement of a qubit. In distinction from the result of a projective measurement, the 
detector outcome is not discrete, since the detector output (for instance, voltage or current) 
is a continuous variable. The outcomes do not even have to correlate with the state of the 
qubit if the detector is uncoupled. Further, the detector variables are subject to noise not 
related to the qubit. Owing to the feedback of the detector at the qubit, this noise affects 
the qubit too. 

In comparison with the text-book projective measurement that instantly provides a result 
and projects the system onto the state corresponding to the result, the CWLM takes time 
both to accumulate the information and to distort the qubit. The time r m required to 
obtain a sufficiently accurate measurement result is called "measurement time" and is a 
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characteristic of a CWLM setup. It is not a duration of an individual measurement in 
this setup: the latter may vary. The distortion is due to the inevitable back action of the 
detector and is characterized by the dephasing rate IV It has been shown^^ that for an 
optimized — quantum limited — detector x m r^ = 1/2 while the "measurement time" r m 
greatly exceeds l/2Y d for less optimal detectors. 

In the context of quantum information theory, CWLM may be understood as an inter- 
action of the qubit with infinitely many ancillary qubits representing the detector degrees 
of freedom. Each ancilla is brought to weakly interact with the qubit for a short time and 
is subsequently measured. Owing to the interaction, the quantum state of the ancillae is 
entangled with the state of the qubit. The detector output is proportional to the sum of 
the measurement results of a large set of ancillae. This allows to transfer quantum infor- 
mation from the qubit to the detector without formal projective measurement of the qubit. 
Therefore the peculiarities of the CWLM can be understood in the framework of a projective 
measurement, although a more complicated one involving the detector degrees of freedom. 
The CWLM can be thus seen as a build-up of an entanglement between the qubit and the 
detector. An outcome of an individual CWLM is the detector output accumulated during 
the time interval of a certain duration r^. Any CWLM can be described as a generalized 
quantum measurement, that involves qubit and detector degrees of freedom. 

The outcome randomly varies from measurement to measurement. We argue here that 
studying statistics of the measurement outcomes of a CWLM is the best way to understand 
and characterize such a measurement. This is especially important for the simultaneous 
measurement of non-commuting variables (say, A and B) we concentrate on in this work. In 
this case, the text-book projective measurement can not help to predict the statistics of the 
results: it would depend on the order of measurements of A and B. This property of the 
measurements in non-commuting bases enables most quantum cryptography^ algorithms 
and has been extensively elucidated in Ref.— . 

One can straightforwardly realize in experiment a CWLM of a quantum system where A 
and B are measured simultaneously. If A and B commute, the statistics of the outcomes of 
sufficiently long CWLM corresponds to the predictions of projective measurement scheme 
(see Sec. III). The projective measurement scheme loses its predictive power if A and B do 
not commute. The reason is that the order of measurement of A and B is not determined 
in the course of a continuous measurement. The statistics of CWLM outputs thus can not 
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be straightforwardly conjectured and has to be evaluated from the quantum mechanical 
treatment of the whole system consisting of the qubit and the detectors. 

In a sharp contrast to the case of commuting variables, the most probable outcome of 
a sufficiently long CWLM of non-commuting variables does not depend on the qubit state. 
Therefore it provides no information about the qubit. The information is however hidden 
in the statistics of random outcomes. Recently, the simultaneous acquisition of two non- 
communing observables was investigated in the framework of CWLM-, and the correlation 
of the random output of two detectors was found to be informative. Not only noise, but the 
whole full counting statistics (FCS) of the non-commuting measurements has been recently 
addressed for an example of many spins traversing the detectors^. 

The structure of the article is as follows. We develop the necessary formalism in Sec. II. 
Our approach stems from the FCS theory of electron transfers^ in the extended Keldysh 
formalism-!, which has been recently discussed^ in the context of the quantum measure- 
ment. At first step, we obtain a Feynman- Vernon action to describe the fluctuations of 
the input and output variables of the detector(s). In the relevant limit, the action is local 
in time. So at the second step we reduce the path integral to the solution of a differen- 
tial equation that appears to be a Bloch-Redfield equation augmented with the counting 
field. In Sec. Ill we exemplify the formalism addressing a relatively simple case of quantum 
non- demolition (QND) measurement^. We evaluate the distribution of the outcomes for a 
single detector and understand the statistics of a recently proposed quantum un-demolition 
measurement^. The main results are presented in Sec. IV where we discuss statistics of mea- 
surement of non-commuting variables for the case of three independent detectors measuring 
the three components of the qubit pseudo-spin. We find the statistical correspondence be- 
tween the three outcomes and three wavefunction components after the measurement. The 
correspondence is characterized by a fidelity that generally increases with the magnitude 
of the outcomes reaching the ideal value 1 in the limit of large magnitudes. Since very 
large outcomes are statistically rare and require long waiting times, this result could be of a 
purely theoretical value. To prove the opposite, we have evaluated the fidelity at moderate 
magnitudes of outcomes and measurement durations and we were able to demonstrate 
the fidelity of 0.95 for ~ 7r m . We term this "quantum monitoring". Ideally, the result of 
the quantum monitoring is a pure state of the qubit and three numbers (detector outputs) 
giving the polarization of the state. The same result can be also achieved by preparing the 
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qubit state of the known polarization, for instance, by a projective measurement along a 
certain axis. The difference is that in the case of preparation the polarization axis is known 
to the observer in advance, while in the case of monitoring it is not so: both the three 
numbers and the state emerge from dynamics of the quantum system that encompasses the 
qubit and the detectors. 

We discuss the relation between the quantum monitoring proposed and the quantum 
algorithms that use the results of partial measurements that we summarize in Sec. V. We 
evaluate the detector action in the Appendix A. We prove in the Appendix B that our 
approach correspond to a Lindblad scheme for a system consisting of the detectors and the 
qubit. 

II. METHOD 

We start the outline of the formalism with a simplest setup where a single detector 
measures a single component of the qubit pseudo-spin. In this case, the Hamiltonian reads 
as follows: 

H — H q + H int + H d , (la) 

3 

H q = ^Hi&i; H int = a 3 Q. (lb) 

i=l 

Here, H q is the Hamiltonian of the qubit generally given by a linear combination of three 
Pauli matrices 6~i (i — 1, 2, 3) corresponding to three components of the qubit pseudo-spin. 
Hi n t gives the coupling between the detector and the third component of the pseudo-spin 
of the qubit, Q being the detector input variable. H d stands for the Hamiltonian of the 
detector. Since we assume linear dynamics of the detector variables, a general form of this 
Hamiltonian is that of a boson bath, 

H d = ^hw k b\b k . 

k 

This encompasses infinitely many boson degrees of freedom labeled by k, b k being the corre- 
sponding annihilation operators. The output of the detector is given by the output variable 
V. An arbitrary linear dynamics is reproduced if both variables Q and V are linear combi- 
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nations of the boson creation/annihilation operators, 

Q = E + , (2) 

it 

V = E + • ( 3 ) 

it 

This is in the spirit of Caldeira-Leggett approach^. In contrast to the work— we do not 
assume thermal equilibrium in the boson bath. In fact, this assumption would be wrong 
for most practical detectors since a signal amplification can not take place in the state of 
thermal equilibrium. The only requirement we impose is that Wick's theorem holds for the 
boson operators involved. This guarantees the linear dynamics of the detector variables. 
Besides, this conveniently allows us not to specify the coefficients Qk, Vk- All information 
about the coefficients and the non-equilibrium boson distribution is incorporated into the 
two-point correlators of the variables explicitly given below (Eqs. 7-8). By virtue of Wick's 
theorem, the averages of all possible products of the detector variables can be expressed in 
terms of these two-point correlators. 

We are interested in the statistics of the detector output variable V. We note that this 
variable is distinct from those of the qubit and in principle even does not have to correlate 
with the qubit state (e.g. if Q — 0). However, since the detector is supposed to measure the 
qubit, there must be a (high) degree of correspondence between the detector output and the 
qubit state. This sets the goal of our calculation: to access the joint statistics of V and the 
qubit variables. 

To achieve the goal, we introduce a counting field x(t) coupled to the output variable V 
and use a modified Feynman- Vernon scheme^ where the evolution of the "bra" and "ket" 
wave function is governed by different Hamiltonians H~ and H + : H ± = H ± Hx{t)V/2. 
± corresponds to two branches of closed time contour respectively^ 1 ^. This scheme was 
first employed in the work—. The counting field x(t) plays a role of the parameter in the 
generating function of the probability distribution of the detector outcomes V(t). This 
generating function is given by: 

Z({ X (t)}) = Tr(Te^f dtH ~R(Q)TeTJ dtH+ ) . (4) 

Tr(- • • ) implying the trace over both detector and qubit variables. Here, T ( T) denotes time 
(reversed) ordering in evolution exponents and R(0) = Pd(0) (^)p(O) is the initial density 
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matrix of whole system. It separates into pd{0) and p(0), the initial density matrix of the 
detector and qubit respectively. This implies that the detector and the qubit do not interact 
before the initial time moment t — 0. 

Next, we employ the path integral representation for the probability-generating 
function^. According to Feynman and Vernon^ (see also^): 

Z({x(t)}) = JvX + VX-e A *Z Iq (Q-,Q + ), (5) 

here X ± (t) are two-dimensional vectors of the detector variables X ± (t) = (Q ± (t),V ± (t)) T , 

VX ± = Y\dQ ± (t)dV ± (t). Q ± (t),V ± (t) are the path-integral variables. Z Iq is called the 

t 

"influence functional" and is given by 

Z Iq (Q-,Q+) = Tr q (Te-H^(^+^Q + W) p (o) 

where Tr q means the trace over qubit space. The action in Eq. (jSj) is bilinear in X^ and 
X to conform to linear dynamics of the detector and will be specified below. The advantage 
of this representation is that the dynamics of infinitely many detector degrees of freedom 
have been reduced to the dynamics of only two relevant fields: Q and V. The influence 
functional describes a non-linear response of the qubit on the fields. 

Let us turn to a specific model of linear dynamics of the detector. Following common 
assumptions about CWLM,— >2 we assume instant detector responses and white (frequency- 
independent) noises. Under these assumptions, a detector is characterized by seven inde- 
pendent parameters: four response functions and three noises. It is convenient to use an 
index i taking values 1 and 2 for input and output variables respectively. With this index, 
we present 4 response functions as a single 2x2 matrix. The noises Sij form a sim- 
ilar matrix. By virtue of Kubo formula, the response functions are expressed in terms of 
expectation values of the operator commutators 



~<[<?(0,Q(O]) 


= a n S(t - t' - 0+) , 


(7a) 


{([Q(tiv(t')]) 


= a 12 5{t - t' - 0+) , 


(7b) 


- h ([V{t\Q{t')\) 


= a 21 5(t -t'- 0+) , 


(7c) 


\([v{tiv{t')\) 


= a 22 5(t - t' - 0+) , 


(7d) 
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where an infinitesimal small positive number + in the <5-function represents small but finite 
response time. 

The noises correspond to the expectation values of symmetrized operator products, 

(ga) 

/(^(0 + ^(«) )) = Wt _ 0> (8b) 

{{ mm±xmt) )) = Si25{t _ tr (8c) 

(( ^ " = S 2i<Ht-0, (8d) 

here, ((AB)) = ((A - (A))(B - (B))) for any operators A and B. 

Let us discuss the physical meaning of the parameters involved. S\\ is the noise of the 
input variable responsible for the back action of the detector and decoherence of the qubit; 
5*22 is the output noise that prevents a fast measurement of the detector outcome. The 
cross-term S 12 = S21 presents the correlation of these two noises. The response function a 2 i 
determines the detector gain: The proportionality coefficient between the detector output 
and the third component of the qubit pseudo-spin, (V) = 021(03). Other response functions 
&i2 ; a 22 ; &ii are respectively related to reverse gain, output and input impedances of the 
detector and are not of immediate interest for us. The detector is characterized with the 
dephasing rate Yd = 2Sn/h 2 and the "measurement time" r m = S^/o^i" 1 ^- The Cauchy- 
Schwartz inequality 



<Sii<S22 — 5"i2 > —(^2i ~ 012) 2 



5! 
4 

imposes an important restrictions on the possible values of the parameters^. Following the 
common assumption, we assume that the reverse gain a 12 is much less than the direct gain 
021 : 021 S> a i2- This condition is commonly required from a good amplifier. Under these 
assumptions, T m Yd > 1/2: one cannot measure a qubit without dephasing it. 
The action A d corresponding to the model reads 

A d = J rft[-^ T (t)(a- 1 ) T 5a~ 1 x(t) + iX T (t)aT 1 x(t) + ix T {t)X(t)} , (9) 

the derivation is outlined in the appendix A. Here we switch to the "quantum" x and 
"classical" X variables defined as follows: x = (X+ - X~)/h, X = (X + + X~)/2. x = 



S 



(0, x) T ■ The 2x2 matrices a and 5 are respectively: 

(10a) 



a ll 0-12 
a 21 a 22 



, 5*11 Si? . 

S = 12 . (10b) 

^21 >->22 

It is important for further advance that the action (Q is local in time. This allows for 
reducing the path integral to a differential equation. The procedure is completely similar to 
the standard reduction of the corresponding path integrals to either Schrodinger or Fokker- 
Planck equations^. One slices time axis into intervals (tk, tk+i) (tk = kAt. k is an integer), 
and takes the path integral in (JSJ) without tracing over the qubit indexes slice by slice. The 
result of the integration at tk is a matrix in qubit indexes, p(tk)- Integrating over x, X in 
the next slice, one finds a linear relation between pt k+1 and p tk : 

p( X ;t k+1 ) = J J] Vx(t)VX(t)e A *S + p( X ;t k )S„; (11) 



tfc<t<tfc+i 



fc+i 



S T = exp ^ ^H q At + a 3 J dtQT(t)j j . (12) 
Since the slice is thin, the exponents may be expanded, 

(tk+l \ tfe+1 

H q At + a 3 j dbQF{t) \ j dtdt'QT(t)QT( t >) + ... 

t k / t k 

and the integration is reduced to evaluation of the averages and the correlators of the fields 
with the action Ad. Collecting terms oc At and taking the limit At — > 0, we obtain a 
differential equation for p(t) that resembles a familiar Bloch-Redfield equation but essentially 
depends on the counting field x^- The resulting equation for p(x; t) reads: 
dp i [T r ~i X 2 (t) „ A ia 2 ix(t),^ . . ^ 

Oj: = -p H q>P\ y~ 22p + — 2 — ( pa3 + a3P > 

—^X{t){ph ~ ozp) ~^-{p- hph) • (13) 

How to apply the equation? Let us consider a single measurement first. Let be the 
duration of the measurement. We collect the detector output during the time interval (0, Td) 
and normalize it by 

V = - [ Td V(t)dt. (14) 
T d Jo 
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To get the statistics of V a , we should assume that x is a constant in the interval (0,7^). 
Indeed, expanding of the generating function in terms of \ gives the averages of products 
of V Q . Let us suppose that the initial density matrix of the qubit is p(0). We solve the 
Eq. (fT51) with the initial conditions p(x ? 0) = p(0). The output is a %-dependent matrix 
p(x) = p(x,t — T d) after the measurement. 

We stress that p(x) appearing in the equation is not the reduced qubit density matrix. 
It is a more interesting and complicated quantity that reflects the joint probability distribu- 
tion of the qubit pseudo-spin components after the measurement and the detector outcome 
collected. To see this, let us define the reduced density matrix of the qubit and the outcome 
V , R(V Q , Vq). It is a matrix in qubit indices and in the outcome values V D , V' Q (see Appendix 
B for the details). Its diagonal elements give the statistics in question. The reduced qubit 
density matrix (with no regard for the value of the outcome) is given by 




the probability distribution of the outcomes (with no regard for the qubit state) reads 

P(V ) = Tr q R(V ,V ), 

and the joint statistics is expressed in terms of the qubit density matrix conditioned to a 
certain value of the output, 

p(V ) = R(V ,V )/P(V ). 

The quantity in use, p(x)? is related to the diagonal elements of thus introduced density 
matrix R by means of Fourier transform^, 

R(V ,V ) = ±_ J dxp(x)e- ixVat (15) 

presenting a generating function for the quasi-distribution R(V , V Q ). Comparing this with 
the above definitions, we find convenient relations 

pit) = p( X = 0), P(V ) = t f ^e-^Tr g p( X ), 

* (v) Jdxe-^Kx) 

It is the main technical advantage of our work that Eq. [13]is similar in form to an elementary 
Bloch-Redfield equation for the qubit density matrix and not much complicated than that 
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one. However, in its augmented form it solves a much more challenging task of finding the 
joint probability distribution of the detector outcome and the qubit state. 

It is important to note that Eq. (ITS]) , as well as Eq. (1171) for multi-detector setup, 
complies with a Lindblad scheme^ 1 ^. We will show this explicitly in the appendix B. This 
guarantees the positivity of the "big" density matrix R(V Q , V Q ). 

The locality in time is a relevant but strong assumption which in fact corresponds to a 
classical detector (indeed, the action (Q does not contain any h) . This is why we do not have 
to worry about possible quantum uncertainties of the detector output that could complicate 
the interpretation of the statistics 28 . 

The scheme can be easily extended to many repetitive (that is, being constantly repeated) 
measurements to comply with the concept of CWLM. Let us consider (infinitely) many 
subsequent measurements. For i-th measurement, the detector output is collected during 
the time interval (t i} t i+ i). This gives a series of outcomes vj . To describe the joint 
statistics, one solves Eq. ( Tl3l) with a piece- wise constant xif)i x(t) = Xi i n the interval 
(ti,ti + i). The solution of the equation at the time moment tM+i depends on M counting 
fields: p = p(xi, ■ ■ ■ ,Xm)- The Fourier transform with respect to all Xi defines the qubit 

(i) 

density matrix conditioned on the outcomes Vo of all M preceding measurements. We 
illustrate two subsequent measurements in the next Section. 

Importantly, the scheme described can also be easily extended to more qubits and/or 
detectors: One just adds extra (counting) fields for detectors and extra Pauli matrices 
for qubits. The case of interest for us is the simultaneous CWLM of three pseudo-spin 
projections of a qubit. The coupling term becomes 

Him = + cr 2 Q 2 + or 3 Q 3 . (16) 

Qk {k = 1, 2, 3) being the input fields of the three detectors. Three counting fields Xk are 
coupled to the corresponding output variables of the three detectors. 

While it is straightforward to write down the equation for general situation, we employ 
a specific model at this point. Namely, we assume for simplicity that the detectors are 
identical and independent. "Identical" implies that the noises and response functions of 
all three detectors are the same. "Independent" implies that no response function relates 
inputs/outputs of two different detectors, neither the noises correlate. Each detector is 
described by the action in corresponding variables. Under these assumptions, the setup is 
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conveniently SU(2) covariant. 

The resulting ^-augmented Bloch-Redfield equation reads 

k=l k=l 

s ^ s ^ 

+ -^^2xk[vk,p] --^(3p-^CT fc po-fc). (17) 

k=l k=l 

Comparing this with Eq. [T3], we see that each detector contributes a term to the equation. 
Each term comes with the corresponding counting field and a-matrix. 



III. STATISTICS OF QND MEASUREMENT 

Before turning to the measurements of non- commutative variables, let us first illustrate 
the formalism with a single-detector setup. Only a single component 03 will be measured. We 
will be interested in a quantum QND setup where successive measurements are performed. 
Such QND measurements have been recently realized for superconducting qubits 35 . To sat- 
isfy the non-demolishing condition^ 3 ., we set H q = ea 3 in ([Taj) so that H q and H int commute. 
In this case, H q can be canceled by transformation to the rotating frame, p — > e tH q t / h p e - tH i t l h 
and will be disregarded from now on. 

Let us perform two measurements that immediately follow each other. During the first 
measurement of duration t\, the detector output is collected in the time interval (0, t\) so the 
measurement outcome is V\ = Jq 1 dtV{t) jt\. Similarly, for the second measurement V2 = 
/ t * 1+ 2 dtV(t)/t 2 . The statistics of the two outcomes is computed from (fT3]) by setting x{t) 
to a piece-wise constant = Xi(X2) during the first(second) time interval and — 
otherwise. To solve the equation, we parameterize t) as follows: 

l + 0"3 / \ 1 — 0"3 

p(x;t) = — — p+{x;t) + — — p-{x;t) 

+&ipi(x;t) + o- 2 p 2 (x;t) , (18) 

where 1 is the unit 2x2 matrix in the qubit space, p± give diagonal elements of the matrix 
and pi ;2 give the non-diagonal ones. Eqs. ffTSl and (|T3j) yield two pairs of separating 
equations: 

_± = i X a 2 ip+ - —S22P+ , (19a) 
dp- x" 2 

— = -i X a 2lP _ - —S22P- , (19b) 
12 



and 

= 7— P2 - (Id + y 5 22 )pi , (20a) 

-QjT = — ^ — Pi - ( T d + y5 22 )p 2 • (20b) 

We first assume that the initial density matrix of the qubit is diagonal. We solve p±{x\ t) 
with initial conditions p±(0) at t = 0. We note that at the initial time t — 0, p does 
not depend on x or Vo an d thus we do not write them explicitly for the initial conditions. 
Solving the equations and Fourier-transforming the generating function yields a very simple 
probability distribution of both outcomes: 

P(^i,v 2 ) = 2^^— p±(0)e 2 e 2 . (21) 

To keep it simple, we have switched here to the dimensionless durations r = Td/r m , and 
outcomes v = V /a 2 i- The result is in fact classical: it does not depend on the dephasing 
rate. In allows for an elementary interpretation: initially, the qubit appears to be either in 
the state + or — , with probabilities p+(0) and p_(0) respectively. The state persists during 
the measurements. The outcome of each measurement is distributed normally around ±1 
with the standard deviation y/T\2 set by the duration of the measurement. We note that the 
persistence of a state is specific for QND measurements, and, as we see in the next Section, 
does not apply to the CWLM of non-commuting variables. 

In Fig. 1, the solid lines show the distribution of outcome P versus v for two different 
durations Ti—2 (long, lower curve) and T\ = 0.3 (short, upper curve). Two obvious peaks 
located at v — ±1 for the line T\ — 2 are due to the two states of the qubit that can be 
distinguished in the course of the long measurement. For the short measurement t\ = 0.3, 
the detector can not resolve the difference between two eigenstates of the qubit. Thus we 
see a single peak at v — broadened by the noise of the output variable. The dotted lines 
show the probability distribution P of the outcome of the second measurement of the same 
duration under condition that the outcome of the first measurement v\ = — 1. For long 
measurement (r 2 = 2), the probability distribution is concentrated near v — — 1. For the 
short measurement (r 2 = 0.3), the distribution is similar to that of the first measurement, 
its average being close to v — 0. This makes a comprehensive illustration of the fact that the 
sufficiently long measurements are repeatable, that is, the result of the second measurement 
is close to the result of the second one. This is not the case for the short measurement. 
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FIG. 1: (color online) Quantum non-demolition setup: Two successive measurements. In each pair 
of the curves, the solid one gives the distribution of outcome of the first measurement while the 
dashed one gives the distribution for the second measurement provided the first measurement gave 
v\ = —1. Lower (upper) pair of curves corresponds to long, = 2 (short, T\ t 2 = 0.3) duration of 
measurement. The long measurement is thus repeatable, the short one is not. 

To illustrate the quantum aspect, let us set the initial density matrix to correspond 
to a pure state with the wave function that is an equal superposition of the base states 
±: p(0) = &\. We are interested in the average value of the corresponding pseudo-spin 
projection a% after the measurement provided the measurement gives the outcome v. We 
evaluate this average if we know the qubit density matrix p{y) conditioned on the outcome 
v. 



<7i (u) = Tr((Jip(f)). 



(22) 



As discussed in the previous Section, this qubit density matrix is computed from the nor- 



malized Fourier transform of p(x), the latter is obtained by solving Eqs. (I20al) and f!20bj) on 
time interval (0, r), r being duration of the measurement. 

Since Eqs. (120bj) do contain the dephasing rate, the result will depend on actual depahs- 
ing. The answer reads 



cos(C 12 vt) _c t 
-e 2 



(23) 



cosh(f r) 

Here, we introduce dimensionless constants C = 4(S , n5 , 2 2 — 5'i2)/(^ a 2i) 2 ~ 1 and C12 = 
2S\2/ha2i- C > characterizes the quality of the detector, C = for a quantum-limited 
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one. C\2 characterizes the correlations of the noises, for a quantum-limited detector Cxi = 
as well. 

Generally, o~i(v, t) quickly decays with increasing r. This manifests the dephasing of the 
superposition by the measurement. Remarkably enough, for a quantum-limited detector 
(C = C\2 = 0) and for a special value of the measurement outcome v = the dephasing is 
absent. The wave function retains its initial value for this particular value of output. This 
fact has been noted in the work— and termed "quantum un-demolition measurement". Let 
us note that the phase shift between the states ±, acquired from the detector, 2 J" 1 " dtQ(t)/H, 
is zero at this (rather improbable^) value of the outcome. We stress that the strict corre- 
spondence between the phase shift and outcome does not hold for a general detector, so that 
ai(v = 0, r) = exp(— Ct/2) decreases with the measurement duration r indicating the non- 
vanishing dephasing of the superposition. We plot <Ti(v; t = 1) versus the detector outcome 
v for the duration of measurement r = 1 in Fig. 2. The solid line is for the quantum-limited 
detector {C = Cyi = 0), and the dotted line is for the a worse detector {C = Cyi — 1). 
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FIG. 2: (color online) "Un-demolition" measurement. The average value of pseudo-spin component 
o"i (v; t) conditioned on the detector outcome v characterizes the dephasing of the superposition after 
a QND measurement of duration r(r = 1 for the plots). A quantum-limited detector (C = C12 = 0, 
upper curve) allows for the quantum "un-demolition" measurement {a\ = 1 at v = 0). This does 
not work for a worse detector (C = C12 = 1, lower curve). 
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IV. THE STATISTICS OF THE CWLM OF NON-COMMUTING VARIABLES 



In the previous Section we provided simple examples to prove the use of the statistical 
approach. Thus encouraged, we turn to the statistics of the CWLM of non-commuting 
variables. The interaction Hamiltonian is now given by Eq. (1161) . We do not want to deal 
with the qubit Hamiltonian H q = ea^ and shall assume that it is removed by transforming 
to the rotating frame. The same transform makes a li2 to rotate with angular velocity e/H. 
To compensate for this, let us presume that the signal from o- 1)2 is collected at frequencies 
e/h rather than at zero frequency as the signal from a 3 is. Mathematically, we define the 
outcomes of the detectors 1, 2 as 

r d Vx, 2 = J " dt[cos(|)V 1)2 (*) t sin(|)F 2il (t)] . (24) 

This can be practically realized in a very same way as it is done in a radio set. One has 
to mix a high frequency signal with a reference signal of the same frequency and detect the 
low-frequency component of the product. 

Let us evaluate /S(xi, X2, Xs)- Without the term H q , the Eq. (fl7|) is readily solved in a 
proper basis in the pseudo-spin space. In this basis, one of the Pauli matrices is defined as 
®x = (X1&1 + X2&2 + X3^3)/Xs, while two others, <r M and a u , are chosen to be orthogonal to 
it. We have introduced Xs as follows: Xs = Vxl + x\ + X%- We parameterize p as follows: 

p = p l + p x a x + p fl (?n + PvVv (25) 

From Eqs. ( 1TTI) and ( l25i) we obtain two equations involving p and p x : 

dpo ■ xl c / oc \ 

-QjT = ^iXsPx ~ , (26a) 

^ = ia 2lX sPo - (2T d + ^-S 22 ) Px . (26b) 

We stress that the CWLM we are about to describe is hardly a measurement of the initial 
state of the qubit. In contrast to QND where the dephasing is limited to 1, 2 components, the 
input variables 0i,2,3 of the detectors randomly rotate the pseudo-spin in all three directions. 
The quantum information about initial state is lost rather quickly: at the time scale of 1/IV 
That is, it is lost before a statistically reliable measurement result can be accumulated. This 
motivates us to choose the unpolarized density matrix as initial condition for the state of 
the qubit before the measurement, 

P(0) = ^1- (27) 
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We will see however that despite the memory lost the CWLM of non-commutative variables 
can be rather informative. 




FIG. 3: (color online) Logarithm of the outcome distribution. The curves from the top to the 
bottom: quantum-limited detector (C<j = 1/2), worse detector (C^ = 2), detector not connected to 
the qubit. 



Let us first discuss the distribution of the detector outputs. We first solve Eqs. (I26a[) and 
(I26bp with the initial condition (12?]) . and then recall Eq. (lol). In the limit of long durations 
r 3> 1, the log of the generating function reads: 



CI 



2 I j^j 

Xs IF 



(28) 



-logZ = r IC d - 

where x« h &s been made dimensionless Xi^22/ a 2i — ► Xi as to give the cumulants of dimen- 
sionless outputs Vi. Here, Cd = r^r m = (C + 1 + C\^)j1 > 1/2. 

The cumulants of the outcomes can be evaluated by taking the derivatives of the log Z 
with respect to Xi- The presence of the qubit enhances output noises of each detector by 
the factor 1 + 1/Ca- 



((ViVj)) 



d 2 log Z . 



r C d 



(29) 



d Xi d Xj ,x " 3 ~ 

There is no correlation of noises between different detectors. Such correlation arise for fourth 
cumulants 

<9 4 logZ. 1 + 2$ 



'I.I 



dXidXj 



2 \XiJ-- 



(C d T) 



(30) 



The distribution is isotropic in three outputs depending on v = yu^ + t> \ + t>f only and in 
the limit of long durations we can calculate it by the saddle-point method, determining an 
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optimal x* corresponding to a given outcome v. We obtain 



\ogP(v) = \ogZ( X *(v)) ] 



d\ogZ 

dxs 



Xs=X* 



= IVT . 



(31) 



X* is purely imaginary. We plot log(P)/r vs. v for three different values of in Fig. 3. 
The solid line is for the quantum-limited detector (C^ = 1/2), the dashed line is for the 
worse detector (Cd = 2). The dotted line is for the detectors not connected to the qubit 
(C<2 = oo). So it is a parabola corresponding to the Gaussian distribution of outcomes in 
this case. We see that the distribution is concentrated at zero. Typical values of outcomes 
v ~ 1/ y/r and for these typical values the distribution can be approximated estimated by a 
Gaussian one P ~ exp (— t> 2 /2(l + 1/Cd))- At larger (and thus atypical) values of outcomes 
(v ~ 1) the distribution is essentially non-Gaussian. We see from the plots that the presence 
of the qubit exponentially enhances probabilities of such outcomes. 

Let us discuss the correlation of the detector outputs and the pseudo-spin after such 
measurement thus turning to the joint statistics of the measurement outcomes and the 
resulting qubit state. 

We characterize the correlation with a fidelity f(v), inner product of the normalized 
vector of the outcomes and averaged pseudo-spin at given outcome v: 



where p(v) is defined in Eq. f fl5|) The fidelity is 1 if the normalized values of the outputs 
precisely give all three pseudo-spin components. Analyzing the saddle-point solution for the 
p(v;r), we obtain that / does not depend on r in the limit r ^> 1. Importantly, at large 
values of the outcomes v 1 the fidelity reaches the ideal value / ~ 1 — C^/v. This, quite 
unexpectedly, enables an efficient quantum monitoring of non-commuting variables. 

The monitoring procedure is as follows. Starting from some initial state, one performs a 
series of repetitive measurements of duration r. The three outcomes of each measurement 
are written down. For most measurements, the values of the outcomes are typical, that is, 
do not exceed the results of such measurements correspond to low fidelity and are therefore 
discarded. One specifically waits for a measurement that gives sufficiently big values of 
outcomes. To decide if the outcomes are sufficiently big, one estimates the fidelity of each 
measurement given the values of outcomes and the relation f(v). If one wants to achieve 
the desired accuracy a^s, one thus waits for the outcomes satisfying f(v) > 1 — a des . The 




(32) 
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big values of outcomes guarantee that the fidelity f(v) is sufficiently high. Sooner or later, 
a measurement gives the sufficiently big outcomes. The quantum monitoring takes place. 
At this moment, the state of the qubit is known with the accuracy desired and it is given 
by the values of outcomes p = diVi/v. Since p 2 = 1, this is a pure state. 

We stress that the monitoring does not constitute a single-shot measurement of all com- 
ponents of the initial unknown quantum state. This would be forbidden by the basic laws of 
quantum mechanics. Indeed, the time required for an accurate monitoring exceeds by far the 
measurement time of the detectors. By this time, the initial state is completely forgotten. 
However, the monitoring gives the observer complete information about the final quantum 
state. 

One also could argue that any pure state of the qubit can be obtained in a simpler fashion. 
One would just choose a proper H q and wait for dissipation to bring the qubit to the state 
of lowest energy. One could also try a projective measurement in a certain basis: after 
several tries, such measurement would give the state desired. We note however that for both 
approaches the resulting pure state is a priori known to the observer. This is not the case 
of quantum monitoring: here we let the quantum system to make "its own choice" of the 
final pure state and do not enforce this choice by any means. 

1 
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FIG. 4: (color online) Fidelity of quantum monitoring / vs. the measurement duration r for 
quantum-limited detector. From upper to lower curve v are respectively 4, 3.5, 3, 2.5, 2, and 1.5. 
Bullets at each curve indicate the value of r at which the probability to get the outcome larger 
than the corresponding v for each curve is 10%; triangles at each curve indicate the value of r at 
which the probability to get the outcome larger than the corresponding v for each curve is 50%. 
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The better the accuracy desired ad es = 1 — / "C 1, the bigger outputs are required, 
v > Cd/ddes- The typical waiting time grows exponentially. To estimate it, we assume the 
duration of each measurement r is in the order of r m . This is because the longer durations 
are not favorable due to decoherence. The waiting time is inversely proportional to the 
probability to have sufficiently high outcomes t w ~ t/P(v; t). Since the success probability 
is exponentially small: P(v;t) ~ exp(— v 2 t/2) from the saddle point solution, we then 
estimate r w ~ exp(t> 2 ). Therefore, we shall expect 

log(r w ) ~ a~ d l s . 

This estimation of the waiting time of successful quantum monitoring sounds pessimistic 
or at least causes a doubt concerning the practical feasibility of the monitoring. To prove 
that the monitoring is practical, we are going to show that a reasonably high fidelity can be 
achieved in a reasonably short time. 

The above arguments are based on the analytical saddle point solution valid at large 
t. Now we investigate the measurements of moderate duration, r ~ 1. This can only 
be done by numerical calculation. Solving Eqs. (I26al) and f!26bj) . and making the Fourier 
transformation according to Eq. (Tl5|) . we evaluate and plot f(v) of a single measurement 
of duration r versus r for the quantum-limited detector and the worse detector {Cd = 2) 
in Fig. 4 and Fig. 5 respectively. Each curve gives /(r) at a given outcome v. At each 
curve, we put a "bullet" to indicate the duration of measurement r at which the probability 
to obtain the outcomes larger than the given one is 10%. Similarly, we put a "triangle" to 
indicate r at which the probability is 50%. To find the positions of the symbols, for each v 
we solve numerically for the values of r that satisfy 

J v + °° Trp(v', t)v' 2 cLv> 



10% (50%) (33) 

J Trp(v', r)v' 2 dv' 

respectively. We see from the plots that / = 0.95 is achieved for a quantum-limited detector 
at v — 4 and r = 0.7. At these parameters, 10% of the measurements are successful, i.e. give 
the output v > 4. To achieve a success, the measurement has to be repeated typically 10 
times given the success probability of 10%. We conclude that the 5% accuracy is typically 
achieved in a time interval ~ 10 x 0.7r m = 7r m . This is not much slower than the QND 
measurement of the same accuracy. 

Let us explain the relation between the monitoring discussed and the quantum algorithms 
that use the results of partial measurement. These algorithms have been introduced in the 
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context of two-qubit systems^^ but also may be applied to a general quantum stated 
Speaking very generally, such algorithms start with a quantum state, pure or mixed, and 
aim at producing another state (pure and/or highly entangled). They proceed in steps. Each 
step involves interaction with ancilla qubits that results in an entanglement of the qubit and 
the ancillae. Importantly, the projective measurement of ancilla qubit (s) is performed at 
each step. The result of this measurement is used to decide upon next step: possible decisions 
include to re-request the initial state, to stop since the fidelity desired is reached, to apply 
a certain quantum gate. 

The quantum monitoring proceeds similarly. It starts with an almost isotropic initial 
state: p = 1. The qubit is being measured by our three independent quantum limit detec- 
tors during a time interval r. The outcome of the detectors is used to make a decision. If 
the sum of the square of the three detectors output V{ of are small: a/v f + v \ + v f < 4, the 
measurement is disregarded: this is an analog of re-requesting the initial state. The measure- 
ment is repeated until the values of the outputs are sufficiently high v = \Jv \ + v f + t>f > 4. 
In this case, the monitoring may stop since its goal is reached: the qubit is purified to a 
state given by the detector outputs v%, vi and v^. Since the Vi (i = 1, 2, 3) are random, the 
purified state of the qubit is also random. Thus one could say that the CWLM monitoring 
implements a quantum algorithm of the sort described, in the same sense as an analogous 
computer may implement discrete computer algorithms. 

V. SUMMARY 

In conclusion, we have shown how to evaluate the full statistics of the outcomes of a 
CWLM on a qubit. We are also able to evaluate joint probability distribution of the outcomes 
and the qubit variables after the measurement. For a single detector, we have illustrated the 
QND measurements and understood the recent proposal of quantum un-demolition measure- 
ment. Most interesting results concern the simultaneous CWLM of three non-commuting 
variables by three detectors. Such "measurement" is obtrusive and typically scrambles the 
initial qubit wavefunction. However, we have demonstrated a high degree of correspondence 
between the wave function after the measurement and the outcomes of the three detectors. 
Therefore, such CWLM may be used for high-fidelity quantum monitoring of the qubit. 
The monitoring in fact amounts to a purification of the qubit state in a random direction 
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FIG. 5: (color online) Fidelity of quantum monitoring / versus the measurement duration r for 
worse detector (C<f = 2). From upper to lower curve v are respectively 4, 3.5, 3, 2.5, 2, and 1.5. 
Bullets at each curve indicate the value of r at which the probability to get the outcome larger 
than the corresponding v for each curve is 10%; triangles at each curve indicate the value of t at 
which the probability to get the outcome larger than the corresponding v for each curve is 50%. 

v = (vi,V2,v 3 )/\v\ at Bloch sphere, t>i_ 3 being random outputs of the detectors. We have 
drawn analogy with quantum algorithms that use the outputs of ancilla measurements to 
decide on the purification degree reached and the quantum gates to be applied. 

The interpretation we give to the results is of course not the only possible one. The com- 
munications with several colleagues have convinced us that "interesting" and "important" 
defy an unambiguous definition as far as theory of quantum measurement is concerned. In 
any case, we have developed calculational tools to access the joint probability of the qubit 
degrees of freedom and the outcomes of linear detectors measuring the qubit. We have also 
derived the representative results. It is up to the reader to conclude. 

Acknowledgments 

H.W. acknowledges the financial support in the framework of NanoNed initiative (project 
DSC. 7023). Y.N. appreciates the participation in 2006 Aspen Summer Program where he 
got the impetus to this work. 




22 



APPENDIX A: DERIVATION OF THE DETECTOR ACTION 



The derivation of the path-integral representation for a set of variables linear in boson 
creation/annihilation is a straightforward task. It is instrumental in dissipative quantum 
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mechanics and therefore is to be found in basic literature on the subject. References 
are usually cited in this respect. Owing to the simplicity of the problem, there are many 
other derivations of the kind that are tailored to specific models (e.g.—) and usually assume 
thermal equilibrium of the boson bath. To avoid any confusion, we present this part of the 
derivation here. We do the derivation in the most general terms possible and specify to the 
concrete model in use at the later stage of the calculation. 

Let Xj be a set of the variables linear in boson creation/annihilation operators, Xj(t) 
being the Heisenberg time-dependent operators of these variables. Let us first disregard the 
coupling with the qubit, so that the time dependence of Xj(t) is governed by the detector 
Hamiltonian only. Explicitly, the Heisenberg equation reads 

-&--h [Hd ' XM - 

We are interested in the generating function of the variables which we present in the 
following form (c. f. Eq. H]) 

Z{{ Xj {t)}) = Tr(Te^^^Wpd(0) Te^ . 

We assume summation over repeating indexes j and skip time indexes of X(t), for 
brevity. Differentiating with respect to the parameters xf(t) of the generating function, one 
reproduces all possible products of the operators Xj(t). 

Let us introduce a path integral over the variables Xp(t) associated with the operators 
Xj(t) by means of the following identity: 



Z({Xj(t)}) = J T>X + VX-Tv 



Tj]5(Xr(t)-X i (t)) e f/^7+ x /)»p d (0) 



t,j 



(Al) 



Here we insert 5-functions that replace the operators Xj(t) by the fluctuating fields X^(t), 
separately for two parts of the Keldysh contour. 
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The use of this representation is that we can treat the coupling between the detector 
and an arbitrary quantum system in the form of the influence functional. If the coupling 
between the detector and the system has the form 

Hint = Xj(t)Sj] 

j 

Sj being operators defined in the subspace of the quantum system, we can formally substitute 
Xj{t) X^(t). The resulting influence functional thus reads 

Z Iq ({X.(t)},{X+(t)}) = Tr(Te-i f dtx i § ip s (0)TeTJ dtx i § >) , (A2) 

where the trace is over the subspace of the quantum system, ps{0) is its initial density matrix 
and the time dependence of Sj is governed by the separate Hamiltonian on the subspace of 
the quantum system. 

Let us turn to the evaluation of the representation (IA1I) . To facilitate the operations with 
5-f unctions, we represent them by means of extra integration over the auxiliary variables 
fc ± (t)//i. At each time moment, 



8{Xf - Xi) = j 



2vr 

With these extra variables, the integral becomes 

Z{{ Xj {t)}) = J VX + VX-Vk + Vk-e^ dtx ^ k ^ + ^ /2) e^ dtx t^t + ^l^ 

Tr [t 'e~ l ■> " dtXjk J p d (0)T ' e~ % ^ ~ dtx J k t] . (A3) 

Let us now take the trace over the boson degrees of freedom. To do this, we use the 
widely known relation 

(eV> = e <^ + ^> 

that holds for A, B that are linear in boson operators under condition of Wick's theorem. 
This allows us to express the trace in terms of the two-point correlators of boson variables, 
(Xi(t)Xj(t')). We may assume (Xj) = without compromising generality. The resulting 
expression reads 

Z({ Xj (t)}) = J VX+VX-Vk+Vk-e^^^^^e^^l^+^^e^- 
A d = - / dtdt' 



^+(t)^ + (0(T(X 4 (t)X,(t'))) + \kr{t)kj{t'){T{Ut)X {t'))) 



+k-(t)k+(t')(X l (t)X J (t')) . (A4) 
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It is instructive to introduce at this stage "classical" and "quantum" variables defining 
X- = Xj ± hxj/2, = Kj ± kj/2h. In these variables, the two-point correlators are 
naturally collected to symmetrized noises t') and Kubo-like response functions Aij(t, t'), 

Sij&t') = ^{X^XjW+XjWXiit)); Aifat) = -|e(t-tO<[*i(t), *;(*')]>■ 
and the action becomes 

A d = - J dtdt' [2i^(0##)-%(M') + iki^Kj^Aij^t')} . 

and does not contain h. One can now integrate over the auxiliary variables Kj(t), kj(t) to 
get the action in terms of Xj (t) , Xj (t) . This Gaussian integral can be readily taken by the 
saddle-point method. For our model, it is convenient to make the time-local approximation 
first. We replace the kernels SV,, Aij with their time-local expressions Eqs.([7],[ED an d perform 
integration over K, k to arrive at Eq. [9j 

APPENDIX B: AUGMENTED BLOCH-REDFIELD EQUATION AND LIND- 
BLAD FORM 

In this appendix we cast the Eq. (fl3l) to the Lindblad form 33 to illustrate the entanglement 
of the qubit and detector and to prove the positivity of the density matrix encompassing 
the qubit and output variable of the detector. 

To specify this underlying density matrix, it is useful to introduce a quantum variable p 
defined as 

p= [ V(t')dt'. 
Jo 

This variable represents the integrated detector output over the interval (0, t). It differs only 
by a factor t from the outcome V Q of the measurement at the same time interval (0, t) as we 
have defined in the main text. 

We denote by p and \p) the eigenvalue and eigenvector of p respectively: p\p) = p\p). The 
subject of our interest is the time-dependent reduced density matrix in the space |spin) (£) \p), 
p(p,p'; t), where "hat" denotes the matrix in the pseudo-spin space. Since p is related to V D 
upon a factor, this matrix is equivalent to R(V , V^) used in the main text. The quantity 
p(x) i n the augmented Bloch-Redfield equation is related to the diagonal part of this density 
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matrix by Fourier transform 

p(p,p;t)= I d X e-™p( X ,t) . (Bl) 



Making the inverse Fourier transform, we deduce from Eqs. (TT5T) and (IB1I) the equation for 
p(p,p;t): 

dp(p,p;t) % S 22 d 2 p(p,p;t) 

--[H q ,p(p,p;t)\ + 



2 



dt /J ' 2 

Q-21 , dpip,P;t) ^ +5 _ dp(p,p;t) , 
2 dp " " dp 
iS 12 , dp(p,p;t) _ ^ dp(p,P',t) \ 



h dp <9p 
—ff(p(p,P;t)-&3P(p,p;t)a 3 ). (B2) 

This is an evolution equation in partial derivatives. 

Next, we demonstrate that Eq. (1B2[) is indeed of a Lindblad type. This also proves that 
the density matrix satisfying the equation has positive diagonal elements. We work in the 
space |spin) (^) \p)- We can check that the operator 6 = i-^ has the following properties: 

(p\po\p) = -i — , (B3a) 

z a .9p(p,p';t) 

\P\op\p) = % q- • (B3b) 
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The Lindblad form of an evolution equation reads 

f = ~\ [H, P ]+l Lt) + [L v , pit]) , (B4) 

V 

where if is a Hermitian operator, generally not coinciding with the qubit Hamiltonian, and 
L v [y = 1, • • • ) are arbitrary operators. We need to prove that a proper choice of the 
Lindblad operators L v and the Hamiltonian H reproduces Eq. 
To do so, we introduce 2 Lindblad operators as follows: 



Ll = V 4^ * s ' (B5a) 

L 2 = \fS^{6 - -ttt-^3 - T^r^s) , (B5b) 
no 22 £o 22 

here ^SnS 22 — 4Sf 2 — h 2 a 21 > is guaranteed by the Cauchy-Schwartz inequality (see Sec. 
II). The Hermitian operator reads: 

H = H q - ^6a 3 . (B6) 
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where H q is the qubit Hamiltonian. We substitute the operators to Eq. (1B4|) to obtain the 
following: 



dp(p,p';t) 1 \ u ~( '4 



dt h 

] S22 , 2 &pM;t) d 2 p(p,p';t) d 2 p{p,p']t) 
2 ^ dpdp' dp 2 dp' 2 



Q21 , dp(p,p';t) . | dp(p,p';t) ^ | . dp(p,p';t) | . dp{p,p';t) , 
iS 12 , dp(p,p';t) ^ dp(p,p';t) ^ _ . dp(p,p';t) _ . dp(p,p';t) s 

^ ^"i^yi^ ^"i^yi /~lrr\' ' 



h dp " dp' " " dp " dp' 

S 

+-^-(fr 3 p(p,p';t)& a - p(p,p';t)) . (B7) 

We made use of the properties (1B3al) and (1B3bl) . We have not done yet, since the above 
equation is for a two-indexed density matrix p(p,p') and not for its diagonal part. We still 
have to prove that the above equation does not mix the diagonal and non-diagonal elements. 
So that, to relate Eq. (1B2I) with Eq. (1B7j) . we change variables as follows: 

p _|_ p> 



2 



(B8a) 



Thus, 



ft = (B8b) 



9 Id d 



9 1, d d 



dp' 2 y dP s dP d ' 
From Eqs. (ET]), (fB9al) and (lB9bl) we obtain: 

dp(Ps,P*t) i r , o S* &p\P.,P A t) 

at = -H^^ p " Pi '^ + -2" — apj — 

a 21 ,dp(P s ,P d ;t), . dp(P s ,P d ;t), 
~( 9P S ° r3 + CT3 9P S } 
iS 12/ dp(P 3 ,P d ;t), , dp(P S} P d ;t). 



(B9b) 



<t 3 - cr 3 - 



h y dP s 00 dP s 
+-^crzp(Ps,P d ;t)a 3 . (BIO) 

Eqs. ( 1B2j) and (IBlOj) has the same form. Therefore, we have proved that p(p,p;t) satisfies 
the Lindblad equation in the |spin) (^) |p) space matrix form and thus is positive. 
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We can straightforwardly extend the above scheme to Eq. (|TTj) of Sec. II that is valid for 
the case of three detectors. In this case, we introduce 6 Lindblad operators 



'4S11S22 - 4^ - H 2 a 2 21 



msL *" (Blla) 

L 2 i = \ /r S^(o i - ~ TT^-^i) > (Bllb) 



where i = 1,2,3. The Hermitian operator reads: 

* = *,-^Z>, (B12) 



i=l 



where 6, = and pi are the eigenvalues of the operator pi = J^Vi(t')dt' (i = 1,2,3) 
corresponding to each detector. 

To conclude, we have shown that the density matrix p(p,p';t) in the space |spin) (^) \p) 
satisfies the Lindblad equation and thus its diagonal part is positive. The diagonal part of 
the matrix is related to the quantity p(x, t) in the Sec. II by the Fourier transform. We note 
that Lindblad operators Eqs. (1B5a| IB5bl IBlla| IBllbl) include both the degree of qubit and 
the degree of the detector (s). This signals the entanglement of the qubit and detector (s) In 
this sense, we measure both the qubit and detector(s) in the context of FCS theory. 



1 V. B. Braginsky and F. Y. Khalili, Quantum Measurement (Cambridge University Press, Cam- 
bridge, 1992). 

2 A. J. Leggett, Science 307, 871 (2005). 

3 V. Cerletti, W. A. Coish, O. Gywat, D. Loss, Nanotechnology 16, R27 (2005). 

4 G. Wendin and V. Shumeiko, in Handbook of Theoretical and Computational Nanotechnology, 
edited by M. Riech and W. Schommers (ASP, Los Angeles, 2006). 

5 M. A. Nielsen and I. L. Chuang, Quantum Computations and Quantum Information (Cambridge 
University Press, Cambridge, 2000). 

6 A. N. Korotkov, Phys. Rev. B 60, 5737 (1999). 

7 D. V. Averin, in Quantum Noise in Mesoscopic Physics, edited by Yu.V. Nazarov (Kluwer, 
Dordrecht, 2003). 

8 A. A. Clerk, S. M. Girvin and A. D. Stone, Phys. Rev. B 67, 165324 (2003). 

28 



9 A. N. Jordan and M. Biittiker, Phys. Rev. Lett. 95, 220401 (2005). 

10 A. N. Korotkov and A. N. Jordan, Phys. Rev. Lett. 97, 166805 (2006). 

11 D. V. Averin, in "Exploring the quantum/classical frontier: recent advances in macroscopic 
quantum phenomena", Ed. by J. R. Frieman and S. Han, (Nova Publishes, Hauppauge, NY, 
2002), p441. 

12 S. A. Gurvitz, Phys. Rev. B 56, 15215 (1997). 

13 A. B. Zorin Phys. Rev. Lett. 76, 4408 (1996). 

14 S. Pilgram and M. Biittiker, Phys. Rev. Lett. 89, 200401 (2002). 

15 D. V. Averin, K. Rabenstein and V. K. Semenov, Phys. Rev. B 73, 094504 (2006). 

16 A. Fedorov, A. Shnirman, G. Schon and A. Kidiyarova-Shevchenko, Phys. Rev. B 75, 224504 
(2007). 

17 N. Gisin, G. Ribordy, W. Tittel and H. Zbinden, Rev. Mod. Phys. 74, 145 (2002). 

18 P. Busch , T. Heinonen, and P. Lahti, Phys. Rep. 452, 155 (2007). 

19 A. Di Lorenzo, G. Campagnano and Yu. V. Nazarov, Phys. Rev. B 73, 125311 (2006). 

20 L. S. Levitov and G. B. Lesovik, JETP Lett. 55, 555 (1992). L. S. Levitov and G. B. Lesovik, 
JETP Lett. 58, 230 (1993). L. S. Levitov, H. Lee and G. B. Lesovik, J. Math. Phys. (NY.) 37, 
4845 (1996). 

21 Yu. V. Nazarov, Ann. Phys. (Leipzig) 8, Spec. Issue SI-193, 507 (1999). 

22 Y. Makhlin, G. Sch6n and A. Shnirman, Phys. Rev. Lett. 85, 4578, (2000); G. B. Lesovik, F. 
Hassler and G. Blatter, Phys. Rev. Lett. 96, 106801 (2006). 

23 V. B. Braginsky and F. Y. Khalili, Rev. Mod. Phys. 68, 1 (1996). 

24 A. O. Caldeira and A. J. Leggett, Ann. Phys. 149, 374 (1983); A. O. Caldeira and A. J. Leggett, 
Physica A 121, 587 (1983). 

25 R. P. Feynman and F. L. Vernon, Ann. Phys. (NY.) 24, 118, (1963). 

26 L. V. Keldysh, Zh. Eksp. Teor. Fiz. 47, 1515 (1964); [Soviet Phys. JETP 20, 1018 (1965)]. 

27 J. Rammer and H. Smith, Rev. Mod. Phys. 58 323 (1986). 

28 Yu. V. Nazarov, M. Kindermann, Eur. Phys. J. B. 35 413, 2003. 

29 M. Kindermann, Yu. V. Nazarov and C. W. J. Beenakker, Phys. Rev. B 69, 035336, (2004), 
appendix A. 

30 Y. Makhlin, G. Schon and A. Shnirman, Rev. Mod. Phys. 73 357, (2001). 

31 H. Kleinert, Path Integrals in Quantum Mechanics, Statistics, Polymer Physics and Financial 

29 



Markets, World Scientific (Singapore), 3rd edition, 2005. 

A. Romito and Yu.V. Nazarov, Phys. Rev. B 70, 212509 (2004). 

G. Lindblad, Comm. Math. Phys. 48, 119 (1976). 

C. Presilla, R. Onofrio, U. Tambini, Ann. Phys. 248, 95 (1996). 

A. Lupascu, S. Saito, T. Picot, P. C. de Groot, C. J. P. M. Harmans and J. E. Mooij, Nat. 
Phys. 3, 119 (2007). 

C. H. Bennett, G. Brassard, S. Popescu, B. Schumacher, J. A. Smolin and W. K. Wootters, 
Phys. Rev. Lett. 76, 722 (1996). 

D. Deutsch, A. Ekert, R. Jozsa, C. Macchiavello, S. Popescu and A. Sanpera, Phys. Rev. Lett. 
77, 2818 (1996). 

H. Nakazato, T. Takazawa and K. Yuasa, Phys. Rev. Lett. 90, 060401 (2003). 



30 



